Atoms versus photons as carriers of quantum states 
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The problem of the complete transfer of quantum states and entanglement in a four qubits system 
composed of two single-mode cavities and two two-level atoms is investigated. The transfer of single 
and double excitation states is discussed for two different coupling configurations between the qubits. 
In the first, the coupling is mediated by the atoms that simultaneously couple to the cavity modes. 
In the second configuration, each atom resides inside one of the cavities and the coupling between 
the cavities is mediated by the overlapping field modes. A proper choice of basis states allows to 
identify states that could be completely transferred between themselves. Simple expressions are 
derived for the conditions for the complete transfer of quantum states and entanglement. These 
conditions impose severe constraints on the evolution of the system in the form of constants of 
motion. The constrains on the evolution of the system imply that not all states can evolve in time, 
and we find that the evolution of the entire system can be confined into that occurring among two 
states only. Detailed analysis show that in the case where the interaction is mediated by the atoms, 
only symmetric superposition states can be completely and reversibly transferred between the atoms 
and the cavity modes. In the case where the interaction is mediated by the overlapping field modes, 
both symmetric and antisymmetric superposition states can be completely transferred. We also 
show that the system is capable to generate purely photonic noon states, but only if the coupling 
is mediated by the atoms, and demonstrate that the ability to generate the noon states relies on 
perfect transfer of an entanglement from the atoms to the cavity modes. 

PACS numbers: 03.67.Bg, 03.67.Hk, 42.50.Dv, 42.50.Pq 



I. INTRODUCTION 

The ability to control the process of transferring quan- 
tum states between remote systems is essential for quan- 
tum information processing and the development of 
quantum computation [1, 2]. The main challenge is the 
achievement of a quantum interface that would be able 
to transfer quantum states with hight fidelity. Several 
schemes have been proposed and experimentally imple- 
mented where a high fidelity transfer of quantum state 
was achieved based on the creation of a strong coupling 
between the systems [3-6]. In this limit, the time scale of 
the transfer process is much shorter than the time scale 
for dissipation in the system due to a coupling to an ex- 
ternal environment. Over this time scale, coherent and 
reversible transfer of a quantum state can be achieved [7] . 

The transfer of a quantum state corresponds to the 
transfer of correlations contained in the states of one 
system to the states of an anther system. In most 
schemes, linear atomic chains or atomic (spin) lattices 
are considered for the transfer of quantum states [8- 
10]. In these models the transfer is mediated by the di- 
rect dipolc-dipolc coupling between neighbouring atoms, 
which induces the flow of an initial excitation through 
the chain [11]. However, the dipole-dipole coupling may 
create correlations between the atoms so the state would 
effectively be transferred through the correlated system. 



* sbougoufTa@hotmail.com 
t zficck@kacst.edu. sa 



Thus, the final state of the system might not be related 
in a simple way to the initial state. The atoms can be 
found in an entangled state even if the initial state was 
a separable state [12]. 

Another schemes involve a chain of remote cavities 
each containing a single atom or optical lattice, the so- 
called Jaynes-Cummings (JC) cells [13-15]. The cells 
could be independent [16-20], or could be coupled to each 
other through the overlapping field modes [21-26], or con- 
nected by short fibres [27-32]. In these schemes, photons 
constitute as carries for the transfer of quantum states 
and controlled transfer is implemented by an appropri- 
ate choice of the coupling strength of the overlapping field 
modes. The cavity modes could be in a vacuum state or 
in a correlated state [33-37]. 

Alternative to the schemes involving photons are mod- 
els in which atoms are treated as a small reservoir or a 
spin bath through which an initial state could be trans- 
ferred to the field modes [38-40] . The atoms forming the 
reservoir can be independent or may interact with each 
other, the two cases that could give different results for 
the transfer process of a quantum state. 

Using the Schrodinger or master equation approach, 
one can study the transfer process of an initial quan- 
tum state for the three classes of systems described above 
(atomic chains or lattices, coupled JC cells, field modes 
coupled through atomic reservoir). In this paper, we ad- 
dress the question of the complete transfer of quantum 
states and entanglement in four-qubit systems composed 
of two single-mode cavities and two two-level atoms. We 
work in the resonance regime and consider two different 
coupling configurations between the qubits specified by 
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two distinctly different types of the interaction Hamilto- 
nians. In the first, we assume that both atoms simul- 
taneously couple to two independent cavity modes re- 
sulting in an effective system analog of the system com- 
posed of two qubits coupled to a small reservoir. Here, 
the coupling between the field qubits is mediated by the 
atoms. The other configuration corresponds to a system 
of two directly coupled JC cells, where the cells are cou- 
pled through the overlapping mode functions of the cav- 
ity fields. In this case, the coupling between the qubits 
is mediated by photons. In what follows, we assume that 
dissipation effects are not present or can be neglected. 
This restricts our results to the case of a strong cou- 
pling between qubits. In fact this is not an overly restric- 
tive limitation regarding the recent progress in the cav- 
ity QED technology, where strong couplings have been 
achieved [41]. 



We shall investigate what kind of quantum states can 
be completely transferred under the action of the two dif- 
ferent interaction Hamiltonians, and under what condi- 
tions the complete transfer of the states corresponds to 
perfect transfer of entanglement between atoms and the 
field modes. We also determine the nature of the states 
associated with the complete transfer of quantum states 
and entanglement. Of course, there is a large number of 
states to which an initial state could be transferred un- 
der the action of the interaction Hamiltonians. However, 
we find that only few of them can be completely trans- 
ferred between themselves. Moreover, we show that the 
complete transfer of a quantum state does not necessary 
mean perfect transfer of entanglement. 



The paper is organized as follows. In Sec. II, we pro- 
vide detailed description of the four-qubit system and the 
coupling configurations between the qubits. A detailed 
calculation of the probability amplitudes of single exci- 
tation states is presented in Sec. Ill, and in Sec. IV we 
extend these calculations to double excitation states. We 
point out the existence of constants of motion that sig- 
nificantly reduce the number of states that could evolve 
in time. We find that the Hilbert space of the system 
splits into independent subspaces, each composed of two 
states only. In other words, the system behaviours in this 
case as composed of independent two-state subsystems. 
No perfect transfer of an entanglement is achieved even if 
the states arc completely transferred between themselves 
when one of the states contains correlations between the 
atoms and the cavity modes. The perfect transfer of the 
entanglement is achieved between states that are factor- 
izable into a product of the atomic and the field states. 
Finally, we summarize our results in Sec. V. Some de- 
tails of the calculation of the logarithmic negativities of 
the atomic and the cavity-field subsystems are presented 
in the Appendix. 



II. GENERAL FORMALISM 

We consider a four qubit system, two identical cavities 
each composed of a single bosonic mode of frequency u) c , 
and two identical atoms modelled as having two energy 
states, a ground state \gi) and an excited state |ej) sep- 
arated by frequency loq and connected by a transition 
dipole moment /ij (i = A,B). The atoms are repre- 
sented by spin lowering and raising operators a~ and crf , 
whereas the cavity modes are represented by the annihi- 
lation and creation operators (a, b), and (a 1 ,^), respec- 
tively. We use indexes (A, B) to label the atoms and (a, b) 
to label the cavity modes. Each of the qubit pairs, either 
atoms (A, B) or the cavity modes (a, 6), can be used as 
a mediator between the qubits. 

Our objective is to determine which of the qubits, 
atoms (fermions) or photons (bosons), are better carri- 
ers of quantum states and entanglement from one qubit 
pair to the other. For this purpose, we concentrate on 
the transfer of an initial state through the system for two 
configurations of the coupling between the qubits. In the 
first configuration, we use atoms as "couplers" by send- 
ing them through the cavities. In the second configura- 
tion, the atoms are assumed to reside inside the cavities, 
with one atom in each cavity, and the coupling between 
the cavities is mediated by overlapping of the evanescent 
mode functions of the cavity modes. Alternatively, this 
could be done by connecting the cavities by a short opti- 
cal fibre. In this kind of coupling, photons constitute as 
carries of an excitation between the cavities. 



(a) 




(b) 
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FIG. 1. (Color online) Schematic diagram of two coupling 
configurations between qubits in the system composed of two 
single-mode cavities (a, b) and two two-level atoms (o\4,o\b). 
In (a), both atoms simultaneously couple to the cavity modes 
resulting in an effective close loop system. In (fe), two Jaynes- 
Cummings cells are coupled with a strength k through the 
overlapping of the mode functions of the cavity fields. 



Figure 1 shows two possible coupling configurations 
in the system. In the first, graphed in Fig. 1(a), each 
atom couples to both of the cavity modes that are not 
directly coupled to each other. In geometrical terms, this 
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configuration is equivalent to a square shape or a closed 
chain configuration of the four qubits. In the second, 
graphed in Fig. 1(b), each atom couples to only one of 
the modes that are directly coupled to each other. The 
coupling between the cavity modes could be arranged 
through the overlap of the mode functions of the cavity 
fields in the region between the cavities. We see from the 
figure that the system is equivalent to an open chain of 
four qubits. This system can be regarded as a double JC 
systems, two JC systems coupled to each other with the 
coupling strength k. 

The two configurations of the coupling between the 
cavities can be distinguished by different forms of the in- 
teraction Hamiltonian of the atoms and the cavity modes. 
The total Hamiltonian of the system can be written as 



where 



H = H a + Hi, (1) 



H Q = fw (aA + &B) + hu) c (tfa + Pb) (2) 



is the free (unperturbed) Hamiltonian of the atoms and 
the cavity modes, and Hi is the interaction between the 
atoms and the cavity modes. The interaction depends on 
the configuration of the coupling between the atoms and 
the cavities and between the cavity modes themselves. 

In the case (a), the interaction Hamiltonian, under the 
rotating-wavc approximation, can be written as 



a+e tkRAB )a 



Hi =h [g a (<t 

+9b (^+a+c lkRAB )b + U.c. 



(3) 



The Hamiltonian consists of two terms. The first term 
describes the interaction of the atoms with the field of 
the cavity a, and the second term describe the interaction 
with the field of the cavity b. The strength of the inter- 
actions is determined by the coupling coefficients g a and 
gb, respectively. In practice we may encounter a situation 
in which atoms travelling through the cavities would be 
at different positions. Therefore, we have introduced the 
phase factor exp(ikRAB) in the coupling coefficients re- 
flecting a phase difference in the oscillation of the atomic 
dipole moments separated by a finite distance Rab- 
If the atoms are separated by less than a wavelength, 
kRAB 1, the phase difference may be taken equal to 
zero that the atomic dipolcs oscillate in phase. When 
the atoms are separated by a half of the wavelength the 
phase difference between the atomic dipoles equals to ir, 
the phase factor exp(ikRAB) = — 1 5 and then the atomic 
dipoles oscillate with opposite phases. 

In the case (&), the interaction Hamiltonian may be 
written as 



Hi = h ygacrxa + g b o B l 
+ Tin (ab f + 5a 1 " J , 



H.c 



(4) 



where (i = a, b) is the coupling coefficient between the 
ith atom and its local cavity field, and k is the coupling 



coefficient between the cavity modes. In the absence of 
the coupling, k = 0, the system described by the Hamil- 
tonian (4) is completely equivalent to two independent 
JC models. 

The calculations for both configurations are based on 
the solution of the Schrodingcr equation for the wave 
function of the system 



ih^- t |*„(t)) = H\v n (t)) 



72=1,2, 



(5) 



where n labels the number of excitations present in the 
system. We will work in the interaction representation 
to remove the time dependence at the optical frequen- 
cies and expand the wave function in the basis of the 
eigenstates of the free Hamiltonian Hq. 

For the case of one excitation (n = 1) present in the 
system, the wave function can be expanded in terms of 
four product state vectors 



l*i(*)>=E C i(*) IJ>' 



(6) 



where | j) are one excitation eigenstates of the free Hamil- 
tonian 

|1) = \e A ,gB,0a,0b) , |2) = |gA,e B ,0 a ,0 6 ) , 

1 3) = \gA,9B,U,0b) , 1 4) = \gA,9B,0 a> h) , (7) 

and Cj (t) is a slowly varying part of the probability am- 
plitude of the state j. The space of the \j) states can 
be formally divided into two subspaces; one involving 
states |1) , 1 2) corresponding to one quantum present in 
either of the atoms, and the other involving states |3) , |4) 
corresponding to the atoms residing in their ground 
states and the excitation present in either of the cavity 
modes or in a superposition of the modes. 

When two quanta of excitation (n = 2) are present in 
the system, the wave function can be expanded in terms 
of eight product state vectors 



|* 2 (t)> = X>(*) li) 



(8) 



where \j) are two-quanta eigenstates of the free Hamil- 
tonian 

|1) = \e A ,e B ,0 a ,0b) , |2) = \e Al gB,Q a , h) , 

1 3) = \gA,e B ,la,0b) , 1 4) = \9A,9B,la,U) i 

1 5) = \eA,gB,U,0b) , 1 6} = \9a,9b, 2 , 0&) , 

1 7) = \g A ,e B ,0a,h), \~8) = \g A ,9B,0a,2 b ). (9) 

The space of the product states can be formally divided 
into two subspaces. One involves states |2) , |3) , |5) , |7) 
corresponding to the two quanta of excitation evenly re- 
distributed over the atomic system and the cavity modes. 
The other set involves states |1) , |4) , |6) , |8) correspond- 
ing to both quanta being either in the atomic system or 
in the cavity modes. 
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III. TRANSFER OF SINGLE EXCITATION 
STATES 

Let us first consider the transfer of single excitation 
states between atoms and the cavity modes. We are 
particularly interested in which states can be maximally 
transferred and whether the complete transfer is accom- 
panied by complete transfer of entanglement. The ques- 
tion of the entanglement transfer will be addressed by 
considering the concurrence of the atomic and the cav- 
ity states. This allows to estimate the efficiency of the 
transfer of the initial state between the atoms and the 
cavity modes that is mediated by atoms, and compare 
the transfer efficiency with that mediated by photons. 



A. Transfer mediated by atoms 

Consider the configuration (a), in which the interac- 
tion is determined by the Hamiltonian (3). By substi- 
tuting the wave function (6) into Eq. (5), we obtain the 
following system of coupled differential equations for the 
probability amplitudes 

iCi = g a C 3 + gbCi, iC 2 = g a C 3 + gbC^, 
iC 3 = -AC 3 +g a (c 1 +C 2 ) , 

iC 4 = -AC 4 +g b (c 1 + C 2 ), (10) 

where A = ujq — uj c is the detuning of the atomic res- 
onance frequency ujq from the cavity frequency u c and 
C 2 = e*C 2 , with e = exp(ifci?AB)- In order to simplify 
the notation without loss of generality, we can assume 
that the coupling strengths g a and g b are both real. 

In order to solve Eqs. (10), we find convenient to in- 
troduce linear combinations of the amplitudes 

W(t) = i[Ci(t) + C 2 (t)}, U{t) = ±=[Ci(t) - C 2 {t)], 

X (t) = ^C 3 (t) + ^C 4 (t), Y(t) = 9 -^C A {t) - ^C 3 (t), 
,9o go go go 

(11) 



where go = \J g\ + g\- The amplitudes (11) correspond 
to orthonormal states 

H = -7= (l e A,.9s) +e\g A: e B )) ® |0 o ,0 6 ) , 
v 2 

\u) = -j= (\eA,9B) - e\g A ,e B }) ® |0 o ,0 6 ) , 

\x) = —{g a \l a ,0 b ) +gb\0 a ,W) ® \9a, 9b), 
9o 

\V) = — (9a\0a,lb) - 9b\la,0b)) ®\9A,gB) ■ (12) 

.9o 

These states are product states of purely atomic and cav- 
ity fields (photonic) states. In the states \w) and \u), the 



atoms are themselves in maximally entangled states, but 
the field modes are in a factorized state. On the other 
hand, in the states \x) and \y), the cavity modes are in 
non-maximally entangled states and the atoms are in fac- 
torized states. This shows that purely atomic entangled 
states can be transferred or converted into purely pho- 
tonic N = 1 noon states [7, 42-45]. 

In terms of the superposition amplitudes, Eqs. (10) 
transform into a simple set of equation 



W(t) = -V2ig X(t), X(t) 
Y(t) = iAY(t), U{t) = 0, 



iAX(t) 



V2ig W(t), 
(13) 



which shows that the amplitudes X(t) and W(t) are 
coupled only among themselves, whereas Y(t) and U(t) 
evolve independently of the others. 

The solutions of Eqs. (13) are readily found to be 



W(t) = W(0)e 
- iX(0) 



hiAt 



A 

cos(r2<) + i— sin(rit) 



-e 2 



X(t) =X(0)c^ lAt 
-iW(0) 



A 



cos(ftt) + i— sin(f2t) 



^■ 90 e^ At sin(ffi), 



Y(t) = Y(0)e 



iAt 



U(t) = £7(0), 



(14) 



where = \j2g^ + A 2 /4 is a detuned one-photon Rabi 
frequency. 

The solutions (14) exhibit a number of interesting fea- 
tures. Firstly, the direct coupling between the ampli- 
tudes Wit) and X(t) indicates that the states \w) and |x) 
form a subspace in the Hilbert space composed of two 
states that can be reversibly transferred between them- 
selves. The remaining amplitudes Y(t) and U(t) each 
evolves independently of the other states, so that the 
states \y) and \u) cannot be accessed from and transferred 
into other states. Once the system is initially prepared in 
either \w) or \x) state, the state vector of the system will 
evolve only between these two states. The system will 
behaviour in this case as a two-state system. Secondly, 
the time evolutions of the atomic and the field states are 
completely symmetric with W(t) and X(t) interchanged, 
despite the fact that the atoms are not equally coupled to 
the cavity modes. Thirdly, the linear superposition U(t) 
does not evolve in time. Thus, if initially U(0) ^ 0, it 
will remain unchanged for all time. In other words, the 
initial population will be trapped in this state and never 
evolve. Fourthly, the time evolution of the superposition 
amplitudes is independent of e, i.e. is independent of the 
relative phase between the atoms. 

The presence of the trapping states has a dramatic 
effect on transfer of entanglement. In order to see it, we 
evaluate concurrences [46] of the atoms and the cavity 
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modes, and find 

Cab® = \{\W{t)\ 2 - \U(0)\ 2 )+2lm[U(0)W*®]\ , 
C ab ®=2\g a9b (\X®\ 2 ~\Y(0)\ 2 ) 

+g 2 a X®Y*® - g 2 X*®Y®\/g 2 . (15) 

These expressions show explicitly that only initial states 
with amplitude U(0) = can be completely transferred 
between the atoms and the cavity modes. However, the 
condition of U(Q) = is necessary but not sufficient for 
the complete transfer of entanglement. It is easy to see, 
since the factor 2g a gb/g^, appearing in the expression 
for C a b(t) is smaller than one, it can be easily verified 
that C a b{tn) could be equal to Cas(O) at some particular 
times t n only if the atoms and the cavity modes have 
the same frequencies (A = 0) and the atoms couple the 
the cavities with the same coupling strengths (g a = gt,). 
Needless to say, U(0) must be zero, and the qubits must 
be indistinguishable for the transfer of the entanglement 
to be complete. Under these conditions an initial, not 
necessary maximally entangled state of the atoms, can be 
completely or perfectly transferred to the cavity modes. 

It is interesting that despite of the complete symmetry 
between W® and X(i), the transfer of maximally en- 
tangled states between the atoms and the cavity modes 
is not symmetric. For example, when the cavity modes 
were initially in the vacuum state but the atoms were 
prepared in the entangled state, i.e. the initial state of 
the system was 

1*0) = ~^f e A9B) + \g A e B )) ® |0 o , 6 ) , (16) 

then the concurrence of the cavity modes varies in time as 

C ab (t)= - l 9 " 9 " sin 2 (to). (17) 

On the other hand, when initially atoms were in their 
ground states but the cavity modes were prepared in the 
N = 1 noon state 

|*o} =^(|l„0 6 ) + |O a l 6 ))<g)|<M,<7B>, (18) 

then the concurrence of the atoms varies as 

Cab® = ff" + fl ff sm 3 (flt). (19) 

z \9a 1 9b) 

Evidently, the variations of the concurrences are not 
equal, and Cab® > C a b® for all times indicating that 
the transfer of the initial maximally entangled state from 
the field modes to the atoms occurs with a better effi- 
ciency than the transfer of entanglement from the atoms 
to the modes. 

The asymmetry in the entanglement transfer is more 
dramatic if the cavity modes were initially prepared in a 
state \x) , (^"(0) = 1). In this case, the concurrences are 

Cab® = sin 2 (fit), C ab ® = - ^f* cos 2 (fit). (20) 

\9a + 9b I 



We see that the atoms, entering the cavities in their 
ground states can be maximally entangled independent 
of whether the cavity modes were initially prepared in 
the maximally entangled state or not. 

The reason for this feature lies in the fact that the sys- 
tem of two atoms unequally coupled to two cavities ef- 
fectively behaviours as a system composed of two atoms 
equally coupled to a single superposition mode. To show 
this, we introduce the collective symmetric atomic oper- 
ator 

»+ = (a+ + ea+) , (21) 

together with symmetric and antisymmetric superposi- 
tion operators of the field modes 

di = (g a a + gbb)/go, d 2 = (gj> - g b aj /g , (22) 

and find that in terms of the superposition operators the 
Hamiltonian (3) takes a simple form 

Hi = hgoafdi + H.c. (23) 

We see that the atoms effectively interact only with the 
symmetric mode d\ . The antisymmetric mode d 2 is com- 
pletely decoupled from the atoms. Thus, in terms of 
the superposition operators, the system effectively be- 
haviours as a single JC system composed of two atoms 
equally coupled to a single mode d\. In terms of the 
state transfer, it means that only symmetric states can 
be transferred between the atoms and the cavity modes. 

B. Transfer mediated by photons 

We now consider the configuration (6) in which each 
atom is located at a fixed position inside one of the cavi- 
ties and the coupling between the cavities is mediated by 
overlapping cavity modes. In this case, an excitation is 
carried by photons. 

Before considering the process of transferring quantum 
states between the qubits, let us first analyze the Hamil- 
tonian of the system in which the interaction is described 
by Eq. (4). Introducing symmetric and antisymmetric 
linear combinations of the atomic and field operators 

&t = {9a°A + 9bvt)/go, 0-+ = {g a ^A - 9bO-s)/9o, 
dt = (fit + gt^ fy/2, a 1 = (V - S f ) /V2, (24) 

we find that the total Hamiltonian of the system can be 
written as 

H = hwo (a \ + crfj) + h(uj c + ^cftd + h(uj c — k)c'c 
+ ^(^-+cta-+H.c). (25) 

The first line of the Hamiltonian (25) represents the free 
energy of the atoms and the superposition modes with 
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the energies of the superposition modes separated by 2k; 
the coupling of the field modes thus lifts their degeneracy. 
The second line of Eq. (25) represents the interaction of 
the superposition field modes with the collective atomic 
systems. 

We note that the structure of the Hamiltonian (25) dif- 
fers significantly from the Hamiltonian we encountered in 
Eq. (23), the counterpart for the coupling mediated by 
the atoms. The most interesting difference is that two 
JC systems coupled to each other by overlapping cav- 
ity modes can be viewed as two independent and non- 
degenerate JC systems, one composed of the symmet- 
ric modes and the other composed of the antisymmetric 
modes. This suggests that in this case both symmetric 
and antisymmetric states could be maximally transferred 
between the atoms and the cavity modes. The transfers 
could occur at two different frequencies, the symmet- 
ric states could be maximally transferred at frequency 
Lu a = uj c + k, whereas the antisymmetric states could be 
transferred at w a = w c - k. 

Since in general the superposition modes are not max- 
imally entangled, maximally entangled states of the field 
modes can be created in the system even if the atoms 
are separable. This feature is opposite of what we en- 
countered in the case of the transfer mediated by the 
atoms, where maximally entangled states between the 
atoms were created with separable cavity modes. 

We now turn to the evaluation of the state vector of 
the system. With the interaction Hamiltonian (25) the 
Schrodinger equation leads to the following equations of 
motion for the probability amplitudes 



iC\ = g a C 3l iC 2 = g b C 4 , 
iC 3 = -ACs + Kd + gad, 
iCi = -AC 4 + kC 3 + g b C 2 . 



(26) 



Equations (26) can be easily solved and interpreted by 
rewriting them as equations of motion for symmetric and 
antisymmetric superpositions 



C s = (C 3 + C 4 )/V2, C a = (C 3 - G0/V2, 

C+ = (g a Ci + gbC 2 )/go, C_ = (g b d - g a d)/go- 



(27) 



They are 



iC s = -(A-k)C s + ^C+, 
V2 

111 II 

iC a = -(A + K)C a + -j=C- + -j=C + , 

,rv - ^Lr + — a 

C+ ~ V2 S+ V2 ai 



(28) 



where u = {g 2 a - g b )/g and w = 2g a g b /g . It is clear 
that the symmetric and antisymmetric modes oscillate 
at different frequencies. Moreover, in the case of identi- 
cal JC systems (u = 0), the modes evolve independently 



from each other. However, the modes may evolve inde- 
pendently even iiu =/= 0. It happens when the field modes 
are well separated in frequency, i.e. when k 3> go- In this 
case, the coupling which exists in general between the 
symmetric and antisymmetric modes is effectively quite 
weak and can be ignored. It is easy to see. By choosing 
a new rotating frame with 

C s = Cse-^-^, C a = C a e^ A+ ^, 

C+ = C+c- 1 ^-^, C- = C-e~ l{A+K)t , (29) 

the equations (28) can than be written in the form 



iC s = ^=C+, 



flip 

V2 

~ w ~ 

iC a = "F-C— 

V2 



V2 



iC+ = (A — k)C+ 



90 -_C, 



V2 V2 



C a e 



2int 



it- = (A + k)C- + ^=c a . 

v2 



(30) 



Clearly, all the u dependent terms are accompanied by 
the exponential factors exp(±2iK<) which in the limit 
of k > go rapidly oscillate in time and thus can be ig- 
nored. This has the effect of decoupling the equations for 
the pair of amplitudes (C S ,C + ) from the pair (C a ,C-). 

Assuming that k ^> go, i.e. that the modes oscillate 
at significantly different frequencies, we then readily find 
analytical expressions for the concurrences, which are of 
the form 



c AB (t) = ^\\c + (t)\ 2 



C ab (t) = \\C s (t)\ 2 -\C a (t)f 



(31) 



where 
C s (t) 



C 2 



\i(A-K)t 



C,(0) cos-fU 



.(A-k)C 5 (0)- V2g C+{0) . 1 

1 K sm 2 QJ 



C+{t) = e ¥(^-*)t 



C+(0) cos^^i 



. (A-K)C + (0) + V2g Q C s (0) . 1 
-1 — sm — \l K t 

n K 2 



(32) 



with Cl K = \j2g~ 1 + (A — k) 2 , and the amplitudes C a (t) 
and C- (t) are obtained from the above by changing 
s — > a, + — > — , k — > —ac, and g — > w. Note that the 
time evolutions of the field states, determined by C s (t), 
and the atomic states, determined by C+(£), are not com- 
pletely symmetric. The evolutions are completely sym- 
metric only at A = «. 

Having the explicit forms of the concurrences we now 
proceed to consider the process of transferring the ini- 
tial maximally entangled states between the atoms and 
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the cavity modes. Provided the atoms are prepared in 
the maximally entangled state (16), the initial conditions 
are C,(0) = C„(0) = 0, C+(0) = (#„ + g b )/(y/2go), and 
C_(0) = (<? a — gb)/(v / 2ffo)- Under this condition, the 
concurrence C ab (t) takes the form 



C ab {t) 



(9a + 9b)' 2 



2gl + (A - «)a 

(5a - 5fc) 2 



y2 ff g + (A-«)2t 



2^ + (A + k) 2 



\/ 2<7q + (A + ft) 2 f 



(33) 



It is clear by inspection of Eq. (33) that in general, the 
concurrence is smaller than one for all times t. However, 
the concurrence could attain the maximal value when 
A = k and g a = g b = <?, in which case C ab (t n ) = 1 at 
the particular times gt n = mr/2, n = 1, 3, . . .. Thus, the 
cells must be identical and A = k for the transfer of the 
maximally entangled state from the atoms to the cavity 
modes to be complete. 

If initially the cavity modes are prepared in the maxi- 
mally entangled N = 1 "noon" state (18), which implies 
C + (0) = C_(0) = C o (0) = and C s (0) = 1, the concur- 
rence CUs(i) is 



C AB (t) 



^9ai 



2cM + (A- K r Sm \2^ 92 ° + ( A - K)H )- 

(34) 



We see that similar to the case of transferring the max- 
imally entangled state from the atoms to the cavity 
modes, the transfer of the maximally entangled state 
from the field modes to the atoms can be complete only 
if 9a = 9b and A = k. 



1- 


w / \ / 




0.8 


'In \ ! s '/ 


i\ 


c °- 6 ~ 


'/ 1 1 \ '7 




t ^ 

° 0.4- 


/ 

'/ i 111 M 


' \ l \ 

■ \ 
1 n 
' i\ 


0.2- 


'/ 1 \ •'/ ' 


« »\ 


n - 


1 1 ; ■ V / \ ' 


> H- 
\ \.» 



3 
gt 




The transfer is far from complete when g a =^ g b and/or k 
close but not equal to A. Notice that independent of the 
relation between g a and gb and whether A = k or not, 
the cavity modes become maximally entangled whenever 
the atoms are disentangled. 

In summary of this section, comparing the results of 
transferring quantum states and entanglement in the con- 
figuration (b) to those of the configuration (a), we find 
that their operations share many common features. For 
example, in both configurations the coupling strengths 
of the atoms to the cavity modes must be equal for the 
transfer of maximally entangled states to be complete. 
However, there are some important differences. In the 
configuration (a), only the symmetric states could be 
completely transferred, whereas in the configuration (b) 
the transfer of both symmetric and antisymmetric states 
is possible. Moreover, in the configuration (a) the com- 
plete transfer occurs at the exact resonance between the 
atomic and the cavity frequencies, A = 0, whereas in the 
configuration (b) the optimal transfers occur at nonzero 
detunings A = ±k. 



IV. TRANSFER OF DOUBLE EXCITATION 
STATES 



Let us now examine the problem of transferring double 
excitation states with two quanta present in the system. 
The two quanta could be in the atoms or in the field 
modes, or could be shared between the atoms and the 
field modes. As one can notice from the form of the 
state vector of the system, Eq. (9), the presence of two 
quanta of excitation increases the number of states in 
which the system can be initially prepared and to which 
it could evolve. Following the same procedure as in the 
previous section, we consider separately the transfer of 
initial states mediated by the atoms and by photons. 



FIG. 2. (Color online) The time evolution of the concurrences 
(a) CAs(i) and (b) C a b(t) for g a = gb = 5, A = « (solid line), 
9a = g, gb = 2g, A = k (dashed line), and g a = gb = 5, A = 
5g, ft = 4g (dashed-dotted line). The system was initially in 
the state (18). 

The above considerations are illustrated in Fig. 2, 
which shows the concurrences Cab(1) and C ab {t) as a 
function of time for several different values of the param- 
eters g a ,gi,,A and k. It is seen that at certain times and 
only for g a = g^ and A = k the entanglement is com- 
pletely transferred from the field modes to the atoms. 



A. Transfer mediated by atoms 

In order to discuss the transfer of double excitation 
states, we make use the wave function (8). Since for the 
transfer mediated by atoms, the best result is expected 
to occur when the atomic transition frequencies are on 
resonance with the cavity field frequencies, therefore we 
set A = in what follows in this subsection. When the 
wave function (8) is inserted into Eq. (5) and we use the 
interaction Hamiltonian (3), we obtain the following sys- 
tem of coupled differential equations for the probability 
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amplitudes 



iD 1 

w 3 

iD 4 
iD 5 

il> 6 



g a (D 3 +D 5 )+g b (D 2 +D 7 ), 
g b Di + g a Di + \/2g b D%, 
g a D\ + gbDi + V2g a D 6 , 
g a (D 2 +D 7 )+g b {D 3 +D 5 ), 
QaDx + g b D 4 + V2g a D e , 
V2g a (D 3 + D 5 ), 



iD 7 = g b Di + g a D 4 + V2g b D 8 
ib s =V2g b (D 2 +D r ), 



(35) 



where A = e*A, (i = 1, 3, 7). 



1. Independent Jaynes-C'ummings systems 

Before continuing with the solution of Eq. (35), we 
briefly comment about the difference between the system 
considered here and that of two independent JC systems, 
extensively studied by other authors [17-20]. When each 
cavity contains one atom, and the cavities are not coupled 
to each other, we then have two independent JC systems. 
In this case, the system of equations for the probability 
amplitudes, Eq. (35), breaks up into three uncoupled sets 
of equations. One composed of four coupled equations 

ibi = g a D 3 + g b D 2 , ib 2 = g a D 4 + g b D 1 , 

ib 3 = g a bi + g b D 4 , ibi = g a D 2 + g b b 3 , (36) 

describing the evolution of the states with a single excita- 
tion present in each of the JC systems, and two separate 
sets each composed of two coupled equations 



and 



iD 5 = V2g a D 6 , tD 6 = V2g a D 5 , (37) 



D 7 = V2g b D 8 , tb s = V2g b D 7 , (38) 



describing the evolution of the states with two excitations 
present in either of the JC systems. 

The set of four equations (36) can be further split into 
two independent sets by introducing symmetric and an- 
tisymmetric superpositions 



D sl = (A + AO/A D s2 = (D 2 + D 3 
D al = (Di - Di)/y/2, D a2 = (D 2 - D 3 )/V2, (39) 

for which Eqs. (36) transform into two separate sets each 
composed of two coupled equations 



iD s i = (g a + gb)D s2 , iD s2 = (g a + g b )D s 



(40) 



and 



The amplitudes (39) correspond to a set of orthonormal 
states 



1*1) = 


-j= (\ca, es, a , 


o b ) 


+ \9A,9B, la, 


h)), 


1*2) = 


(\eA,gB,0a, 


h) 


+ \gA,e-B,la, 


Ob)), 


al) = 


{\sa, es, a , 


o b ) 


- 5A,.9B, la, 


h)), 


a2) = 


A= (\eA,9B,0a, 


h) 


- \9A, es, l a , 


Ob))- 



We see that despite of the complexity, the dynamics of 
the system simplifies to that occurring among states of 
four independent doublets 

\sl) o \s2) , \al) o \a2) , |5) O |6) , \7) o |8) . (43) 

This shows explicitly that the Hilbert space of the system 
splits into four independent subspaces each composed of 
two states only. In other words, the system behaviours 
in this case as composed of four independent two-state 
subsystems. The states of different subsystems oscillate 
at different Rabi frequencies. For example, the states of 
the doublet (\sl) , | s2) ) oscillate at the Rabi frequency 
9a + gb- The states of the remaining subsystems oscil- 
late at frequencies g a — gb, V^ga, and \[2g b , respectively. 
Note that in the limit of g a = gb, the states of the dou- 
blet (|al) , |a2)) decouple from each other and become 
trapping states. 

The states (42) are entangled states of the atoms with 
the field modes. However, in these states both, the atoms 
and the field modes arc themselves disentangled. It is 
easy to see. Properties of the atoms alone can be de- 
scribed by a reduced density operator obtained by tracing 
the total density operator p(t) over the cavity modes 

p AB (t) = Tr ab p(t) = (0 a b \ P {t) |0 a b ) 

+ (l o 6 | p(t) \i a o b ) + (0 a l b \ p{t) \0 a l b ) 
+ (l a l b \ p(t) \l a l b ) + (2 a b \ p{t) |2 a D ) 
+ {0 a 2 b \p{t)\0 a 2 b ). (44) 

It is easy to show that in the basis spanned by four state 
vectors, \e A e B ) , \&AgB) , IffAes) , |ffA5s>, the density op- 
erator pab (t) is diagonal with all off-diagonal elements 
(coherences) equal to zero. Consequently, there is no 
possibility for entanglement between the atoms. 

Yonac et al. [17, 18] have proposed to include an auxil- 
iary zero-excitation state, |0) = \gA, gB,0 a ,0 b ) that per- 
mits the two-photon coherence to be involved in the 
atomic dynamics. When the state |0) is included, the 
following density matrix pab (t) is obtained 



PAB(t) = 



iD a i = (g a - g b )D a2 , iD a2 = (g a - gb)D a i- (41) 



pu(t) P u(t) 

P22 {t) 

p 33 (t) 

041 (*) Mi) 



(45) 
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where Pll (t) = \D 1 (t)\ 2 ,p 22 (t) = \D 2 (t)\ 2 + \D 5 {t)\ 2 7 
Paa(t) = \D 3 (t)\ 2 + \D 7 (t)\ 2 , PiA (t) = \D 4 (i)\ 2 + \D & {t)\ 2 + 
\D 8 {t)\ 2 + \D (t)\ 2 , Pu (t) = Di(t)D (t), and £>„(*) is the 
amplitude of the auxiliary state |0). 

The coherence may result in an entanglement between 
the atoms. It is easy to see, the concurrence evaluated 
from the matrix (45) takes the form 



C AB (t) = 2max{0,C 1 (t),C 2 {t)}, 



where 



Ci(f) 
Ca(*) 



-Vpn(t)p4A(t), 
\ P u(t)\-\D 2 (t)\ 2 



\D 3 (t)\' 



(46) 



(47) 



It is seen that C 2 (t) can be positive. Thus, the inclusion 
of |0) into the atomic dynamics is critical for entangle- 
ment between two atoms placed inside independent cav- 
ities. 

We may summarize that the Hilbert space of the in- 
dependent JC systems can be spanned in terms of four 
independent subspaces, each composed of two states that 
can be exchanged between themselves during the evolu- 
tion. In the case of equal coupling strengths of the atoms 
to the field modes, one of the subspaces is composed of 
trapping states that do not evolve. 



2. Coupled Jaynes-Cummings systems 

We now proceed to solve the set of Eqs. (35), which 
shall allow us to discuss the transfer efficiency of doubly 
excited states in the system composed of two coupled JC 
systems. Since the transfer process can be strongly af- 
fected by the presence of trapping states in the system, 
we solve Eqs. (35) by finding constant of motion that 
determine states which do not evolve in time. This will 
allow us to predict which initial states can be completely 
transferred in the system. 

We first observe that D 7 = D 2 and b 3 = D5. This 
suggests the introduction of symmetric and antisymmet- 
ric superpositions of the states with the excitation shared 
between the atoms and the cavity modes 



D f 



2, D s2 = {D 2 + D 7 )/V2 



D a i = (-D3 - D 5 )/V2, D a2 = (D2-D 7 )/V2, (48) 

and find that then the set of Eqs. (35) splits into two 
separate sets, one involving six coupled differential equa- 
tions 



iD x = V2{g a D sl +g b D s2 ), 
iD sl = V2{g b D 1 + g a D 4 + V2g b D 6 ), 
W s2 = V2(g a D 1 + g b D 4 + V2g a D 8 ), 

W4 = \/2(g a D s2 + g b D s i), 

iD 6 = 2g a D sl , 

iD s = 2g b D s2 , 



and the other composed of two constants of motion 

ib al = 0, ib a2 = 0. (50) 

Equations (49) are cumbersome due to the asymmetry 
in general between the coupling constants of the atoms 
to the field modes, g a ^ g b . However, we have seen in the 
case of single-quantum states that the complete transfer 
of the states and the optimal transfer of entanglement 
occur when the atoms couple to both cavities with the 
same strengths, g a = g b . Therefore, in what follows we 
shall restrict to the case of g a = g b = g. 

When g a = g b , additional constants of motion appear. 
To see this, we introduce symmetric and antisymmetric 
superposition of the probability amplitudes 



D u = {D sl + D s2 )/V2, D w = (D s i - D s2 )/V2, 

D p = (D 6 + A0/V2, D q = (D 6 - A0/V2, (51) 

we find that the set of coupled equations (49) can be 
split into two separate sets, one involving four coupled 
equations 

iDi = 2gD u , 

ib u = 2gD 1 + 2gD 4 + 2gD p , 
W 4 = 2gD u , 



iDp = 2gD u , 



(52) 



and other involving two coupled equations 

it) w = 2gD q , iD q = 2gD w . (53) 

The set of equations (52) can be reduced further to a set 
of two coupled equations 

iD u = 2V3gD z , iD z = 2V3gD u , (54) 

and two constants of motion 

w m = 0, ib n = 0, (55) 

where 

D m = {Da - D p )/V2, D n = (D x - D 4 )/V2, (56) 
and 

D z = {bt +D 4 + D p )/ VS. (57) 

According to the above predictions, we may conclude 
that the Hilbert space of the system can be split into 
three independent subspaces, one subspace composed of 
four trapping states 



l al > = -J| i\9A,e B ) 



\eA,gB}) ® [la, Ob) 



\a2) = -j= (\e A ,gB) - \9A,e B )) ® |0 o , h) 
\m) = I (V2\l a , lb) - \2 a , b ) - |0 a , 2 b 



(49) 



\9a,9b) 



\n) = —j= (\e A ,e B ,0a,0 b ) - \g A ,9B, la, lb)) , (58) 
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and two subspaces composed of doublets 
, , 1 



w 



\q) 



(\g A ,e B ) + \e A ,g B )) ® (|lo,0 6 ) - |0„, 1 6 » , 
(|2 a ,0 b )-|0 a ,2 b ))®| ffj4 , 3B ), (59) 



1 



V2 



and 



(|su,e B > + \e A ,9B)) ® (|la,0&) + |0 aj h)) , 



v/3 



[|e^, e B ,0 a ,0 b ) + Ifl'AjS'B, la, It) 



-;^(|2 a ,0 6 ) + |0 a ,2 6 ))®l<M,SB) 



(60) 



The evolution of the doultes (59) and (60) is deter- 
mined by the amplitudes (D Wl D q ) and {D u , D z ) 1 respec- 
tively. Their time evolutions arc readily found to be 

D w {t) = D w (0) cos (2gt) - D q (0) sin (2gt) , 

D q {t) = D q (0) cos (2gt) - D w (0) sin {2gt) , (61) 

and 

D u (t) = D u {0) cos(2V3gt) - iD z (0) sin(^2V35 
D z (t) = D z (0) cosfoVSgt) - iD u (0) sin^VSgt) . (62) 

Thus, due to the presence of four trapping states, the 
evolution of the system confines to that occurring only 
among states \w) <-> \q) and \u) <-> \z). Each pair of 
states evolves in time independently of the other, and 
the pair (\w) , \q)) oscillates at frequency 2g, whereas the 
pair (\u) , \z)) oscillates at frequency 2\/3g. 

Several features of the dynamically evolving states (59) 
and (60) are worth noting. One can see that the state \w) 
is a product state of a single-excitation entangled state 
of the atoms and the asymmetric TV = 1 noon state of 
the cavity modes [7, 42-45]. The state \q) is a product 
state of the atomic state in which both atoms are in their 
ground states and the cavity modes being in the asym- 
metric TV — 2 noon state. This shows that an initial 
state in which atoms are prepared in the maximally en- 
tangled state and the cavity modes are prepared in the 
antisymmetric TV = 1 noon state can be converted to the 
asymmetric TV = 2 noon state. In other words, in this 
scheme we can achieve a fully dynamical generation of a 
purely photonic TV = 2 noon states. 

The other pair of states, \u) and \z), involves the sym- 
metric TV = 1 and TV = 2 noon states. One can notice 
that in contrast to the asymmetric TV = 1 noon state, 
an initial symmetric TV = 1 noon state cannot be com- 
pletely converted to the symmetric TV = 2 noon state. It 
is rather converted into the four-qubit state \z), which 
involves a superposition of the TV = 2 noon state and 
two other states with the excitation evenly redistributed 
between the atoms or between the cavity modes. 



The reason for the difference in the transfer properties 
is in the trapping effect of the single excitation states. Ac- 
cording to Eqs. (12) and (13), the antisymmetric TV = 1 
noon state is a trapping state whereas the symmetric 
state evolves during the atomic transit time through the 
cavities. Thus, the simultaneous evolution of the atomic 
and the symmetric noon states results in a redistribu- 
tion of the excitation between the atoms or between the 
cavities. 

Another interesting feature of the states (59) and (60) 
is that the complete transfer requires the presence of an 
initial entanglement of the cavity modes. In other words, 
the complete transfer of an initial state arises from the 
quantum nature of the field. It is particularly seen in 
the form of the reduced density operator p AB {t), which 
written in the computational basis is not diagonal 



p A B (t) 



Pn(t) 

P2a(*) P2 3 (t) 

psa(*) 

p 4 i(t) 



(63) 



where Pll (t) = ^(t)] 2 , p 22 (t) = \D 2 (t)\ 2 + \D 5 (t)\ 2 , 
(to® = \D 3 (t)\ 2 + \D 7 (t)\ 2 , p 44 (t) = \D 4 (t)\ 2 + \D 6 (t)\ 2 + 
\D s (t)\ 2 , and p 23 (t) = D*(t)D 5 (t) + D 2 (t)D*(t). 

The matrix (63) is of an X-state form due to the pres- 
ence of the coherence p 2 3,{t) between the atoms. If the ini- 
tial state of the system is taken to be one of the states (59) 
or (60), then none of the trapping states (58) is popu- 
lated. As a result, the coherence p 2 z{t) takes a form 

P2s(t) = l(\D u (t)\ 2 + \D w (t)\ 2 ). (64) 

We see that the coherence is completely determined by 
the probabilities (populations) of the states \u) and \w). 
Since the states involve maximally entangled single- 
quantum states of the atoms and the field modes, this 
implies that the coherence can be different from zero only 
if either the atoms or the field modes are in entangled 
state. 

The coherence is necessary for entanglement between 
the atoms. The concurrence evaluated from the ma- 
trix (63) is given by 



C AB (t) = 2max{0,C 1 (t),C 2 (t)}, 



where 



Ci(t) = \ P23 (t)\ - Vpn(iW(t), 
C 2 (t) = -(\D 2 (t)\ 2 + \D 3 (t)\ 2 ). 



(65) 



(66) 



Clearly, C\ (t) can be positive indicating an entanglement 
between the atoms. Thus, an entanglement can be gener- 
ated between the atoms without the need of introducing 
the auxiliary state |0). 

We now present some numerical calculations of en- 
tanglement evolution between the atoms and the cavity 
modes to check if the complete transfer of the double ex- 
citation states is accompanied by the complete transfer 



11 



of entanglement. We evaluate logarithmic negativities 
Nab and N a i, for three different initial states. A detailed 
description of the method of evaluating the logarithmic 
negativities is given in the Appendix. 




FIG. 3. (Color online) Time evolution of the logarithmic neg- 
ativities NAB(t) (blue dashed line) and N a i(t) (red solid line) 
for the initial state \w), Eq. (59). The other parameters are 
A = and g b = g a = 3- 

Figure 3 shows the logarithmic negativities Nab^) 
and N a b(t) as a function of the dimensionless time gt 
for the initial state \w) . According to Eq. (61), at times 
t = mr/(4:g), n = 1,3,..., when the state \w) is com- 
pletely transferred to the state \q), the atoms are ex- 
pected to be separable and the field modes should be 
found in the maximally entangled N = 2 noon state. We 
see that at the initial time t = 0, both negativities are 
maximal, Nab(0) = N a b(0) = 1, indicating that initially 
the atoms and the field modes were maximally entangled. 
As time progresses, the atoms become periodically disen- 
tangled but the field modes never completely disentangle. 
At times t = nir/^Ag), when the atoms are separable, the 
modes again become maximally entangled. However, the 
maximal entanglement of the modes does not mean that 
at these times the modes returned to their initial state. 
At times t = mr/(Ag), the system is in the \q) state, in 
which the atoms arc in their ground states and the field 
modes arc in the asymmetric N = 2 noon state. 

If the initial state is not one of the states (59) or (60), 
only a partial transfer of the initial entanglement may 
occur. To illustrate this, we show in Fig. 4 the time 
evolution of the logarithmic negativities when the system 
starts from an initial state 

|*o) = ^d 1 ^) + |0«0 fc » ® \g A g B ) , (67) 

in which the cavity modes are in the maximally entan- 
gled two-photon state and the atoms are in their ground 
states. We see that in this case the concurrence N a b(t) 
never becomes zero as time develops and NAB(t) never 
reaches the maximal value of NAB(t) = 1. Thus, no com- 



plete transfer of the entanglement occurs. The modes re- 
main highly entangled at times the atomic entanglement 
reaches its maximum value. This result is as expected 
from the above simple discussion that for states different 
than the states (59) and (60), a part of the initial exci- 
tation is trapped in one of the trapping state preventing 
the complete transfer of the initial state to occur. 




FIG. 4. (Color online) Time evolution of the logarithmic neg- 
ativities NAait) (blue dashed line) and N ab (t) (red solid line) 
when the atoms in their ground states |<?a<?b) pass through 
the cavities prepared initially in the maximally entangled 
state (|l a lj,) + |0 a 0f,))/\/2. The other parameters are A = 
and g b = g a = g. 

It is interesting to contrast the initial state (67), which 
involves the auxiliary zero-excitation state, with an ini- 
tial state involving states containing two quanta of ex- 
citation. The former was treated in the case of inde- 
pendent JC cells [17-20] and shown crucial for the cre- 
ation of en entanglement between the atoms. Figure 5 
shows the time evolution of the logarithmic negativities 
for an initial state in which the atoms are in their ground 
states and the field modes arc in the N = 2 noon state 
(|2 Q b ) + |0 a 2 6 ))/v / 2. While Fig. 4 shows that an initial 
entanglement of the field modes can be transferred to the 
atoms when the modes are prepared in the state (67), 
Fig. 5 shows no such transfer occurs when the modes are 
initially prepared in the symmetric N = 2 noon state. 



B. Transfer mediated by photons 

We now turn into the other configuration in which the 
coupling between the cavities is mediated by the over- 
lapped cavity modes. Such coupling configuration, il- 
lustrated in Fig. 1(b), is represented by the interaction 
Hamiltonian (4). What we are going to calculate is essen- 
tially the same as for the first scenario. We use the same 
basis states, Eqs. (9), to find the equations of motion for 
the probability amplitudes required to study the dynam- 
ics of entanglement. As in the previous case, the starting 
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FIG. 5. (Color online) Time evolution of the logarithmic neg- 
ativities NAB(t) (blue dashed line) and N a b(t) (red solid line) 
when two atoms in their ground states pass through the cav- 
ities prepared initially in the symmetric N — 2 noon state. 
The other parameters are A = 0, gb = g a = g. 



point is the Schrodinger equation (5), which with the 
Hamiltonian (4) leads to the following set of equations of 
motion for the probability amplitudes 

ib x = AD 1 +g b D 2 +g a D 3 , 

iD 2 = 9bDi + g a D± + kD 5 , 

«A = g a Di + g b Di + nD 7l 

W 4 = -AD 4 + g a D 2 + g b D 3 + V2k(D 6 + D 8 ), 

W 5 = V2g a D 6 + nD 2l 

W 6 = -AD 6 + V2kD 4 + V2g a D 5 , 

W 7 = kD 3 + V2g b D s , 

it) s = -AD 8 + V2nDi + V2g b D 7 . (68) 

Before proceeding with the solution of the equa- 
tions (68), we first rearrange them to see if there are 
any trapping states in the system. For such a treatment 
it is convenient to introduce linear combinations 



b 2 


= (A" 




^3 


= (D 2 


-A)/ 


V2, 


A 


= (D 5 - 




b 7 


= (D B 


-A)/ 


V2, 


D 6 


- (A;- 


h£» 8 )/V2, 


As 


= (A 


-A)/ 


V2, (69) 



for which Eqs. (68), with g a = g b = g, form two separate 
sets, one involving five coupled equations 

ibx = AD i + V2gD 2 , 
iD 2 = V2gD 1 + V2gD 4 + kD 5 , 
iD A = -A£>4 + V2gD 2 + 2k A, 
ib 5 = V2gD 6 + Kb 2 , 

W 6 = -AD 6 + 2kD 4 + V2.gZ) 5 , (70) 



and the other involving three coupled equations 
W 3 = kD 7 , 
W 7 = nD 3 + V2gD 8 , 
W s = -AD S + V2gD 7 . 



(71) 



It is easily verified that the determinants of the coeffi- 
cients of the coupled differential equations (70) and (71) 
are different from zero only if A ^ 0. Thus, in the case 
of A 7^ 0, there are no trapping states in the system. 



1. The case of A = 

Let us first specialise Eqs. (70) and (71) to the case of 
exact resonance, A = 0, in which trapping states occur. 
As before, the occurrence of trapping states is identified 
by the presence of constants of motion. It is easily verified 
from Eqs. (70) that a linear superposition 



2 L> 4 - V2gnD 5 + (k 2 - g 2 )D } 



(72) 



in which fl 2 = \/2g 4 + k 4 , is a constant of motion, 
i.e. ib a = 0. 

Similarly, the system of coupled equations (71) can be 
rewritten in terms of linear superpositions 



Dr. 



V2~9f, K h 



A = £d> 



V2g~ 



(73) 



where D,' = \/2g 2 + n 2 , for which we find that 
W a = 0, 

iD 7 = SI' Dp, iDp = Sl'D 7 . 



(74) 



Thus, the amplitude D a is a constant of motion whereas 
the amplitudes b 7 and Dp undergo a simple sinusoidal 
evolution 

D 7 (t) = D 7 (Q) cos Sl't- A(0) sin Sl't, 

Dp(t) = Dp(0)cosSl't- £> 7 (0)sinfi't. (75) 

Evidently, states corresponding to the combinations A 
and Dp can be completely transferred between them- 
selves with frequency SI' . 

The states corresponding to the combinations A 
and Dr arc of the form 



1 7) = (\e A ,9BAa,0b 



9A, es, a , l b )) , 
7 (\eA,gB,0 a , U) - \gA,e B ,l a ,0 b }) 



V2SI 

(\2 a ,0 b ) - \0 a ,2 b ))®\g A ,g B ) , 



(76) 



from which one can easily notice that the states contain 
correlations between the atoms and the field modes. This 
implies that the complete transfer of the states will not 
lead to the complete transfer of entanglement. 
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2. The case of A ^ 

Let us now examine the situation when A ^ and 
focus our attention on the case of g a = g b = g. With 
nonzero detuning the dynamics of the system becomes 
more complicated and simple analytical results for the 
probability amplitudes are not possible. However, in the 
limit of the strong coupling between the cavities, k~^> g, 
one can uncover some interesting results. By introducing 
symmetric and antisymmetric combinations 

A, = {Di + D 6 )/V2, D c = {D 4 - D 6 )/y/2, 
D x = (D 2 + D B )/y/2, D d = (D 2 - D 5 )/V2, (77) 
and going into a rotating frame through the relations 

D v = A,c- t(A - 2K)t , Dx = Dxe-^- 2 ^, 
D t = 5 f e-< A+2K > ( , D# = Ae- t(A+2K)t , 
Di = Die^ 8A *, (78) 

we find that Eqs. (70) can be written in terms of these 
combinations as 

i£) n = V2gD x , 

iDx = (A - k)Dx + V2gD n + V^gDxe 21 ^, 
iD t = V2gD^, 

ib<, = (A + n)D. g + V2gD e + V2gb!e- 2lKt , 

iD x = g (D x c 2lKt + D#c- 2iKt ) . (79) 

We see that the terms proportional to b\ play the 
role of a coupling between the symmetric and antisym- 
metric combinations. These are accompanied by terms 
cxp(±2int) that oscillate in time with a frequency 2k. 
In the limit of k > 5, we can make a secular approxi- 
mation, in which we ignore the rapidly oscillating terms 
cxp(±2ifti). This has the effect of decoupling the equa- 
tions of motion for the symmetric combinations 

W v = V2gb x , 

iDx = (A - k)D x + V2gb, v (80) 

from the equations of motion for the antisymmetric com- 
binations 

iD, = V2gD n , 

iDv = (A + k)A> + V2gD e , (81) 

that the amplitudes of the symmetric and antisymmetric 
combinations evolve independently. 

It is seen that the symmetric combinations oscillate 
with unequal frequencies that differ by A — k, whereas 
the antisymmetric combinations oscillate with frequen- 
cies that differ by A + n. The difference of the frequencies 
will lead to oscillations of the amplitudes that are not pe- 
riodic. This imperfection would result in an incomplete 



transfer of the states. The transfer could be complete 
only if the amplitudes oscillate with equal frequency. It 
is seen from Eq. (80) that for A = k, the symmetric com- 
binations oscillate with equal frequencies which results in 
sinusoidal oscillations of the amplitudes 

D v (t) = D v (0) cos{V2gt) - Dx(0)sm(V2gt), 

Dx{t) = b x {0) cos(V2gt) - D n (0) sin(\/2gt). (82) 

Similarly, when A = — n the antisymmetric combinations 
oscillate with equal frequencies and then their time evo- 
lutions are perfectly sinusoidal 

D e (t) = D e {0) cos{x/2gt) - b a (0) sm(\/2gt), 

D#(t) = A?(0) cos(V2gt) - D e {0) sm{\/2gt). (83) 

Clearly, states corresponding to the symmetric combina- 
tions and those corresponding to the antisymmetric com- 
binations of the probability amplitudes can be completely 
transferred between themselves. 

The states corresponding to the symmetric combina- 
tions D v and Dx can be written in the form 

|r?) = - (V2|l a ,l b ) + |2 a ,0 6 ) + |0 a ,26))®|.gA,ffB>, 

I A) = \{\9A,e B ) + \e A ,g B ))®{\la,Ub) + KAb)) > (84) 

whereas the states corresponding to the antisymmetric 
combinations D e and D$ can be written as follows 

|e) = i(\/2|l OJ l 6 ) - |2 a ,0 fa ) - \0 a ,2 b )^®\g A ,g B ) , 

|0) = i(|<u, e B ) - \e A ,g B })®(\la, Ob) - |0 Q , 1 6 » . (85) 

Notice an important property of the states (84) and (85) 
that they do not contain correlations between the atoms 
and the field modes. The states are of the form of product 
states of the atomic states and those of the field states. 
This is in distinct contrast to the states for A = 0, 
which contain correlations between the atoms and the 
field modes. An another interesting property of the states 
(84) and (85) is that the atoms may completely disentan- 
gle during the evolution of the states, contrast to the field 
modes that can never disentangle. This is illustrated in 
Fig. 6 which shows the time evolution of the logarithmic 
negativities when the system was initially prepared in the 
state | A). We see that the initial entanglement between 
the atoms and that between the field modes recurs pe- 
riodically. The atomic entanglement periodically falls to 
zero in contrast to the entanglement between the cavity 
mode that never ceases. 

In closing this section, we would like to point out that 
the two coupling configurations between the cavities lead 
to quite different results for transferring double excita- 
tion states. Comparing Eqs. (59) and (60), the quantum 
states that are completely transferred when the transfer 
is mediated by the atoms with Eqs. (84) and (85) for 
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FIG. 6. (Color online) Time evolution of the logarithmic neg- 
ativities Nab(£) (blue dashed line) and N ab (t) (red solid line) 
for the case in which the system was initially prepared in the 
symmetric state |A). The other parameters are gt — g a = g 
and A = k = lOg. 



the states that arc maximally transferred when the cou- 
pling is mediated by photons, we see that the coupling 
mediated by the atoms may result in a " nonclassical" 
N = 2 noon state \q), whereas no such state is created 
by the coupling mediated by photons. In the latter, two 
states are created, \r]} and |e) that resemble very much 
a two-photon "classical" state created at an asymmetric 
beam-splitter with unequal reflection and transmission 
coefficients [47]. 



CONCLUSIONS 



The problem of the complete transfer of quantum 
states has been examined in a four-qubit system com- 
posed of two single-mode cavities and two two-level 
atoms. We have considered two different coupling con- 
figurations between the qubits specified by two distinctly 
different types of the interaction Hamiltonians. In the 
first, the coupling is mediated by "flying" atoms that si- 
multaneously couple to the field modes of the cavities, 
whereas in the second configuration, each atom resides 
inside one of the cavities and the coupling between the 
cavities is mediated by the overlapping field modes. 

We have demonstrated that the problem of the com- 
plete transfer of states is indeed tractable, that a proper 
choice of the basis states for the state vector of the sys- 
tem allows to confine the dynamics of the system to that 
occurring between two states only. In other words, we 
have found that the time evolution of the state vector of 
the system can be divided into two state pairs that evolve 
in time plus time stationary trapping states. 

Moreover, we have found that the complete transfer of 
a quantum state does not necessary mean perfect trans- 
fer of entanglement. An initial entanglement can be per- 



fectly and reversibly transferred between the atoms and 
the field modes only if the states are in the form of prod- 
uct states of atomic and field states. If the states contain 
correlations between the atoms and the field modes, no 
perfect transfer of the entanglement is possible. 

We have also shown that apart from the similarities 
in the dynamics of the system for the two configurations 
there are some important differences. In particular, the 
complete transfer of quantum states is more manifested 
in the system where the coupling is mediated by the over- 
lapping field modes. In this case not only the symmet- 
ric but also the antisymmetric states can be completely 
transferred between themselves. However, a purely pho- 
tonic N = 2 noon state can be created only if the coupling 
is mediated by the atoms. 



Appendix A: Logarithmic negativities of the atomic 
and field systems 



In this Appendix we provide some details of the eval- 
uation of logarithmic negativities Nab and N ab of the 
atomic and cavity-mode subsystems for the case when 
two quanta of excitation are present in the system. The 
logarithmic negativities are evaluated from the knowl- 
edge of the reduced density matrices of the subsystems. 

The reduced density operator pab of the atomic sub- 
system is obtained by tracing the total density operator p 
over the cavity modes 

Pab = Tr ab (p) = (Q a O b \p \0 a b ) + (l a b \ p\l a O b ) 
+ (0 Q l 6 |p|0 Q l 6 ) + (l a l b \p\l a l b ) 
+ (2 a b \ p \2 a b ) + (0 a 2 b \ p \0 a 2 b ) . (Al) 

In the basis spanned by four state vectors, 
le^es) , \e A g B ) , \g A e B ) , \9a9b), the density opera- 
tor pab is of the form 



PAB 



( Pn 



V o 





P22 PlZ 
P32 P33 



\ 




p A J 



where p u = \D 1 \ 2 ,p 2 2 



\D, 



Pa 



I£sl 2 , 



P23 



i P33 = 

DID, 



(A2) 



^3| ' 

D 2 D* 7 . 



The logarithmic negativity reads 



N AB {t) = max{0, log 2 [l - 2 j ui(t)],log 2 [l - 2p 2 {t)}} , 

(A3) 

where 



M2O) = 5 [Pll (*) +P44(t) - VW (t) -P44(t)] 2 +Mp2 3 (t)\ 2 

(A4) 

Similarly, the reduced density operator for the cavity 
modes is obtained by tracing the total density operator 



15 



over the atoms, and is given by 

Pab = Tr AB (p) = (e A e B \p\e A e B ) + {e. A g B \ p\e A g B ) 
+ (9Ae B \p \gAe B ) + {9a9b\ P \9a9b) ■ (A5) 

The determination of entanglement between the cavity 
modes is a slightly more complicated than for the atoms. 
The reason is that in the case of two excitations present 
in the system, the subspace of the cavity modes is panned 
by six states, |0 a h ) , |0 Q 2 b ) , |l a l b ) , |2 a h ) , |l a 6 ) , |0„l t ). 
In this basis, the reduced density operator p ab has the 
block diagonal form 



Pab 



( Pll 

















\ 





P22 


P23 


P2A 













P32 


P33 


P3A 













P42 


P43 


Pa 
























P56 


/ 


V o 











A35 


P66 



(A6) 



where p xl = \D 1 \ 2 ,p 2 2 = \D 8 \ 2 ,p 3 3 = \D 4 \ 2 ,p 44: = 
W,p 55 = \D 3 \ 2 + \D b \ 2 lPm = \D 2 \ 2 + \D 7 \ 2 ,p 23 = 
Ph = D 8 Dl,p 2 4 = p% 2 = D 8 D£,p 34 = p% 3 = D A D* & and 
P56 = P* 65 = D*D 5 + D 3 D* 7 . 

The logarithmic negativity N a b reads 



iV a6 = max{0,log 2 (l + 2^| Atl |)}, (A7) 



where pi are the negative eigenvalues of the partial trans- 
pose matrix p a b Ts ■ 



M. A. Nielsen and I. L. Chuang, Quantum Computation 

and Quantum Information, (Cambridge University Press, [22 

Cambridge, 2000). 

J. I. Cirac, P. Zoller, H. J. Kimble, and H. Mabuchi, [23 
Phys. Rev. Lett. 78, 3221 (1997). 

C. Monroe, Nature (London) 424, 839 (2003). [24 
T. Wilk, S. C. Webster, A. Kuhn, and G. Rempe, Science 
317, 488 (2007). [25 
A. D. Boozer, A. Boca, R. Miller, T. E. Northup, and H. 
J. Kimble, Phys. Rev. Lett. 98, 193601 (2007). [26 
M. Paternostro, G. Adesso, and S. Campbell, Phys. Rev. 
A 80, 062318 (2009). [27 

G. Nikoghosyan, M. J. Hartmann, and M. B. Plenio, [28 
Phys. Rev. Lett. 108, 123603 (2012). 

I. Kamleitner, J. Cresser, and J. Twamley, Phys. Rev. A [29 
77, 032331 (2008). [30 
L. Vinet and A. Zhedanov, Phys. Rev. A 85, 012323 [31 
(2012). 

Z. M. Wang, R. S. Ma, C. A. Bishop, and Y. J. Gu, Phys. [32 
Rev. A 86, 022330 (2012). [33 
C. J. Mewton and Z. Ficek, Phys. Rev. A 71, 060311(R) [34 
(2005). 

Z. Ficek and R. Tanas, Phys. Rep. 372, 369 (2002). [35 

E. T. Jaynes and F. W. Cummings, Proc. IEEE 51, 89 
(1963). [36 

H. Paul, Ann. der Phys. 466, 411 (1963). [37 
M. Brune, F. Schmidt-Kaler, A. Maali, J. Dreyer, E. Ha- 

gley, J. M. Raimond, and S. Haroche, Phys. Rev. Lett. [38 
76, 1800 (1996). 

I. Sainz and G. Bjork, Phys. Rev. A 76, 042313 (2007). [39 
M. Yonac and J. H. Eberly, J. Phys. B: At. Mol. Opt. 
Phys. 39, S621 (2006). [40 
M. Yonac and J. H. Eberly, J. Phys. B: At. Mol. Opt. 
Phys. 40, S45 (2007). [41 
M. Yonac and J. H. Eberly, Phys. Rev. A 82, 022321 
(2010). [42 
S. Chan, M. D. Reid, and Z. Ficek, J. Phys. B: At. Mol. [43 
Opt. Phys. 42, 065507 (2009); 43, 215505 (2010). 

F. K. Nohama and J. A. Roversi, J. Mod. Opt. 54, 1139 



(2007). 

C. D. Ogden, E. K. Irish, and M. S. Kim, Phys. Rev. 
A78, 063805 (2008). 

F. A. A. El-Orany and M. R. B. Wahiddin, J. Phys. B: 
At. Mol. Opt. Phys. 43, 085502 (2010). 
Y. L. Dong, S. Q. Zhu, W. L. You, Phys. Rev. A 85, 
023833 (2012). 

Y. Guo, J. Li, T. Zhang, and M. Paternostro, 
arXiv: 1206.3477 (2012). 

P. Xue, Z. Ficek, and B. C. Sanders, arXiv: 1204. 1680 
(2012). 

T. Pellizzari, Phys. Rev. Lett. 79, 5242 (1997). 

A. Serafini, S. Mancini, and S. Bose, Phys. Rev. Lett. 96, 
010503 (2006). 

P. Peng and F. L. Li, Phys. Rev. A 75, 062320 (2007). 
S. B. Zheng, Appl. Phys. Lett. 94, 154101 (2009). 
S. B. Zheng, Z. B. Yang, and Y. Xia, Phys. Rev. A 81, 
015804 (2010). 

L. H. Sun and G. X. Li, Phys. Rev. A 85, 065801 (2012). 

S. Bougouffa, Optics Commun. 283, 2989 (2010). 

R. Tahira, M. Ikram, S. Bougouffa, and M. S. Zubairy, 

J. Phys. B: At. Mol. Opt. Phys. 43, 035502 (2010). 

S. Bougouffa and A. Hindi, Phys. Scr. T143, 014006 

(2011) . 

S. Bougouffa, AIP Conf. Proc. 1370, 192 (2011). 

S. Bougouffa and Z. Ficek, Phys. Scr. T147, 014005 

(2012) . 

F. Fassioli Olsen, A. Olaya-Castro, and N. F. Johnson, 
Journal of Physics: Conference Series 84, 012006 (2007). 
C. E. Lopez, G. Romero, and J. C. Retamal, Phys. Rev. 
A 81, 062114 (2010). 

W. Shu and T. Yu, J. Phys. B: At. Mol. Opt. Phys. 44, 
225501 (2011). 

R. J. Thompson, Q. A. Turchette, O. Carnal, and H. J. 
Kimble, Phys. Rev. A 57, 3084 (1998). 

B. C. Sanders, Phys. Rev. A 40, 2417 (1989). 

A. N. Boto, P. Kok, D. S. Abrams, S. L. Braunstein, C. 
P. Williams, and J. P. Dowling, Phys. Rev. Lett. 85, 2733 
(2000). 



16 



[44] I. Afek, O. Ambar, and Y. Silberberg, Science 328, 879 
(2010). 

[45] H. Wang, et al, Phys. Rev. Lett. 106, 060401 (2011). 
[46] K. Wootters, Phys. Rev. Lett. 80, 2248 (1998). 



[47] L. Mandel and E. Wolf, Optical Coherence and Quantum 
Optics (Cambridge: Cambridge University Press, 1995). 



